o 

(N 

< 
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We study the pinning of a two-dimensional membrane to a patterned substrate within elastic 
theory both in the bending rigidity and in the strain dominated regimes. We find that both the 
in-plane strains and the bending rigidity can lead to depinning. We show from energetic arguments 
that the system experiences a first order phase transition between the attached configuration to a 
partially detached one when the relevant parameters of the substrate are varied, and we construct a 
qualitative phase diagram. Our results are confirmed through analytical solutions for some simple 
geometries of the substrate's profile. We apply our model to the case of graphene on top of a 
Si02 substrate and show that typical orders of magnitude for corrugations imply graphene will be 
partially detached from the substrate. 
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I. INTRODUCTION 



Until recently, the study of two dimensional (2D) mem- 
branes was developed mainly for its theoretical interest 
and its applications to biological systems which could be 
well approximated by the 2D membrane model, as well as 
soft matter systems E. Nowadays, however, with the ex- 
perimental discovery of graphene 1^"^ (a two dimensional 
graphite sheet), we have in our hands the opportunity of 
studying a truly 2D membrane. It has been proven that 
the membrane aspect of graphene, and in particular the 
presence or not of a substrate, plays an essential role to 
characterize its behavior ISHSl Graphene presents intrin- 
sic ripples 121, inherent to its 2D nature, which can inter- 
act with the propagating electrons and affect transport 
properties J^. In most experimental settings up to date, 
though, graphene is deposited on top of a substrate, ei- 
ther purposely patterned or presenting random disorder. 
A relevant question then is to determine how the spatial 
structure of the substrate affects that of graphene. This 
kind of study also opens the possibility of controlling the 
properties of graphene by patterning appropriately the 
substrate. Experiments have shown that the morphology 
of a graphene membrane on top of a substrate is largely 
determined by the substrate's profile 1^^^, as opposed to 
suspended graphene. The attachment of graphene to a 
corrugated surface leads to the bending and stretching of 
the graphene layer, so that the depinning of the layer may 
become energetically favorable. For device-construction, 
as well as for the interpretation of experimental data, it 
is important to know if the graphene sheet is completely 
pinned to the substrate or if there are regions for which 
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depinning occurs and graphene is suspended. 

Given its experimental relevance, in this work we ad- 
dress the problem of determining which is the stable con- 
figuration of a membrane on top of a substrate which 
presents either depressions or protrusions. Although we 
will treat the problem in the context of graphene physics, 
our results are general. We analyze this problem from a 
general field theory framework, in which we show the 
possibility of a phase transition between a pinned con- 
figuration to a partially depinned one, where relevant 
parameters of the patterned substrate act as control pa- 
rameters. We turn then to analyze some simple substrate 
geometries which allow for analytical solutions. We will 
show that these examples quantitatively confirm our phe- 
nomenological, qualitative model. Our work is the first 
in the graphene literature that takes into account the ef- 
fect of in-plane strains for detachment of the membrane. 
We show that there is a length scale for the substrate's 
pattern beyond which the in-plane strains are dominant 
and can lead to depinning. This length scale marks the 
crossover from a regime in which the bending rigidity of 
the membrane is dominant energetically. 

In what follows we will analyze the depinning of a 
membrane from the substrate for the two different limit- 
ing regimes mentioned above '221. Firstly, in Section pl| we 
introduce the model for the free energy of a membrane 
on top of a substrate and by means of scaling arguments 
we establish the possibility of a phase transition for the 
system between two possible stable equilibrium configu- 
rations: the membrane being completely attached to the 
substrate, or otherwise it being partially detached. From 
this we are able to construct a qualitative phase diagram 
for the system. In the following sections we proceed to 
a quantitative analysis for a given geometry of the sub- 
strate profile. We consider a substrate with a Gaussian 
depression or protuberance and we obtain analytical so- 
lutions for the two limiting regimes, the bending rigidity 
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FIG. 1: Pictorial representation of a membrane on top of a 
random substrate, partially conforming to the substrate. The 
height of the membrane is represented by a field /i(x) while 
the top surface of the substrate is represented by a field s(x), 
as discussed in the main text. 



dominated regime in Section |III| and the strain domi- 
nated regime in Section |TV] In both cases, we show that 
the system presents a first order phase transition from 
pinned-to-depinned as the ratio of width to height of the 
substrate's profile is varied. A discussion and possible 
experimental consequences are presented in Section IV] 



II. MODEL AND QUALITATIVE PHASE 
DIAGRAM 

We consider a tethered membrane which lies on top of 
a substrate. We use the de Monge parametrizationfiSl^ by 
which the membrane is parametrized by (x, /i(x)), where 
h is the height with respect to some reference plane and 
X = [x, y) are the in-plane coordinates. In the same way, 
the profile of the substrate is represented by (x, s(x)), as 
shown schematically in Figure [l] 

We assume as a first approximation that the membrane 
couples to the substrate through its out-of-plane modes 
(also denominated fiexural modes), via a contact force 
characterized by a surface tension 75. Previous works 
that study the attachment of a membrane to a substrate, 
have used the so called Deryagin approximation, which 
approximates the interaction potential between the mem- 
brane and the substrate as a harmonic potential ^^. How- 
ever this approximation results in a strongly confining 
potential. In our case we are interested in studying the 
stability of the pinned configuration and the possibility 
of detachment, and therefore a contact force approxima- 
tion is more appropriate. Moreover, the interaction be- 
tween graphene and a substrate has been studied in Ref. 
[T7]and it has been shown that the attractive interaction 
force decays as the inverse distance to a power that de- 
pends on the type of interaction (a power of 2 in the 
case of undoped Si02). In that work it was also shown 
that the coupling strength decays roughly four orders of 
magnitude when the graphene sheet is not pinned to the 
substrate. These considerations justify the use of a con- 
tact force which is finite when graphene is conforming 
to the substrate and zero otherwise. This is, of course, 
an idealization of our model since we are disregarding 
the equilibrium distance between the substrate and the 
membrane, which for graphene on a Si02 substrate is of 
the order of 5 A12. The free energy for the membrane on 



top of the substrate within this approximation is given 

by 
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where /i and A are the Lame coefficients and k is the 
bending rigidity of the membrane, and S is the sur- 
face of contact between the membrane and the substrate. 
Throughout this paper we will use the accepted values of 
the elastic and bending parameters for graphene at room 
temperature. The bending rigidity is given by k « 1 eV 
^^, and the Lame coefficients are given by /i « 10 eV 
A-2 and A « 2 eV A'^ 22). We take the value of the 
coupling constant strength as 75 = 2 meVA , corre- 
sponding to the maximum est imate d pinning strength for 
graphene on a Si02 substrate l^liffl. The functional depen- 
dence of F\u, h, s] on the substrate's profile field s(x) is 
given implicitly through the contact term, being /i(x) = 
s(x) when the membrane is attached to the substrate. 
(V^/i(x)) is the local mean curvature of the membrane 
and the local intrinsic curvature is encoded in the strain 
tensor '^^t Uij — ^ {diUj + djUi + dihdjh), with u(x) the 
in-plane phonon modes and the i, j = 1,2 index the two 
components of the field. Since the action is quadratic in 
these modes, they can be integrated out 23 to obtain an 
effective free energy e~'^^ff^'^'^^ = / Pu e""'^'""'''*! with 
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the transverse projector and where 



we have used the expression for the Young modulus in 
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2 , ^ . We are interested in analyzing the 
possible detachment of the graphene sheet from the sub- 
strate, and in particular, to find the configuration which 
is energetically favorable. The general procedure would 
be to minimize the free energy Eq. ([2| given a profile of 
the substrate to find the stable solution. However, the 
non-linearity of Eq. ^ makes this program impossible to 
follow analytically, even for the most simple geometries. 
We are then obliged to make use of approximations if 
we are to make any analytical progress. It is usually as- 
sumed that the in-plane stresses are small and therefore 
their contribution, encapsulated in the quartic order term 
of the effective energy Fi,ff[h, s], can be neglected. How- 
ever this is true only if the height ffuctuations are not too 
big, as we proceed to show. If we consider a substrate of 
average height ffuctuations S over a length scale i, from 
Eq. ([1]) we see that the bending energy of a membrane 
attached to this substrate scales as 
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where we have used that V ^ L~^ and the area J d?x ~ 
L^. On the other hand, by similar arguments (note that 
Uij ^ 5^/L^), the elastic energy due to in-plane strains 



is given roughly by 
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(4) 



Therefore the elastic energy due to in-plane strains is the 
main contribution to the total energy of the membrane 
if EanS" 3> K. This analysis is valid except for quasi one 
dimensional (ID) geometries, where the height profile of 
the substrate is constant along one direction. For this 
case it is easy to show that the in-plane strains are com- 
pletely screened by the height fluctuations and hence the 
in-plane stresses are zero, and the only contribution to 
the elastic energy is due to the bending rigidity. 

With the previous analysis we have then arrived to a 
length scale 



E, 



(5) 



that determines a crossover from a bending rigidity dom- 
inated regime (BD regime) for S < I, to a, strain domi- 
nated regime (SD regime) for S > l. With the values for 
the elastic parameters of graphene given above, C w 1 aS2! 
Note that this scale is of the order of magnitude of the 
lattice spacing and in principle this would imply that the 
BD regime for graphene is greatly suppressed '^■^ . How- 
ever recent atomistic simulations have shown that ther- 
mal height fluctuations of this magnitude are possible I^Zl. 
Moreover, the study presented in Ref. "22" shows that the 
continuum model can still be applied in th is limit. This 
scale can thus be realized in grapheneI2l2l and therefore 
the crossover is of experimental relevance. 

We can now study the two limiting regimes separately. 
For the BD regime, the free energy for the membrane can 
be approximated by: 
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To solve for the equilibrium conflguration, we look for the 
saddle point solutions of equation (|6l with a partially de- 
tached membrane and study their stability. Minimizing 
with respect to the height /i(x) yields the bi-harmonic 
equation within the detached region: 



(v2)';i(x) = o 



(7) 



to be solved together with the appropriate boundary con- 
ditions, while h{x.) = s(x) in the pinned region. The 
boundary conditions have to be imposed at the boundary 
of the surface S, that is the curve at which the membrane 
starts to detach from the substrate. If we parametrize 
this closed curve by x* = dS, the boundary conditions 
are given by: 



V/i(x*) = Vs(x* 



(9) 



The curve x* itself is unknown, and can be determined 
by an extra boundary condition which implies a disconti- 
nuity in the second derivatives due to the surface tension 
force at the curve of detachment x* 123 
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Alternatively, it is equivalent to flnd the extreme of the 
free energy Eq. (l6]) as a function of x*. For a given, ar- 
bitrary proflle of the substrate s(x), it is to be expected 
that the free energy will have many extrema, correspond- 
ing to unstable and metastable conflgurations. The curve 
Xg corresponding to a global minimum will give the sta- 
ble equilibrium conflguration, if this is the null curve then 
the stable conflguration is the totally pinned membrane. 
We see then that the curve x* emerges as a natural order 
parameter of the problem, between two possible states 
of the system: a null curve Xq = corresponding to a 
membrane which is completely attached to the substrate, 
and a flnite value of the function Xq which gives a par- 
tially detached membrane. A scalar order parameter can 
be obtained, for example, by taking the total length of 
the curve |x* | H2l. Our analytical results for the particu- 
lar geometries studied, to be developed in the following 
sections, show that the pinned conflguration is always 
at least a metastable minimum and hence the pinned- 
to-depinned transition is always of flrst order. We can 
argue this has to be true in general for smoothly corru- 
gated substrates as follows. If we consider a small de- 
viation of the system from the totally attached conflgu- 
ration |xq| — 0, described by a small detachment curve 
I Sx* I , the energy cost due to depinning is proportional 
to the minimal area enclosed by the curve, ~ |^x*| . On 
the other hand, the smoothness of the substrate implies 
that, for small enough |(5x*|, the area delimited by this 
curve is locally flat and hence the gain in energy due to 
the relaxation of bending and stretching of the membrane 
is negligible. Hence the pinned conflguration is always a 
local minimum of the energy and the phase transition to 
a partially detached conflguration is of flrst order, due to 
the development of new metastable states with the varia- 
tion of the control parameters. It is safe to assume that, 
for flxed external conditions, these control parameters 
will be related to the characteristic width and height of 
the substrate's corrugations. To simplify the analysis, we 
can consider the problem of a single depression or protu- 
berance in the substrate. Intuitively it is to be expected 
that the stability of the pinned conflguration, given a 
coupling strength 75 and bending rigidity k, will depend 
on the aspect ratio of the substrate's proflle. A simple 
energetic argument gives an estimate for this threshold. 
The interaction energy between the graphene layer and 
the substrate in a region of area L^ is 
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/i(x*) = s(x*) 



(8) 



while, as we saw previously, the bending energy cost of 
height corrugations of scale S is given by Eq. (Isl. The 



change between the regime where the pinning energy is 
dominant and the layer is attached to the substrate, to 
the regime where the cost in bending energy leads to the 
detachment of the layer, is governed by the ratio 
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(12) 



The membrane will prefer to attach to the substrate in 



the limit 



> 1, which translates into a condition for 



the substrate profile: 



L2 



< 



(13) 



indicating that pinning is favored for shallower depres- 
sions. 

Within the bending rigidity approximation, detach- 
ment can occur due to the high bending energy cost that 
competes with the energy gain due to pinning. In the op- 
posite regime, S > I, the in-plane stresses are dominant 
and we should consider the possible detachment due to 
these modes. For this case the free energy (fl| can be 
approximated by: 
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(14) 

In this limit, the approximate free energy given by 
Eq. (14) still contains the non-linear coupling between 



the in-plane and out-of-plane modes and hence further 
approximations are necessary for obtaining analytical re- 
sults '32. We will introduce these approximations in Sec- 
tion [TV] when we solve the system for a particular geome- 
try of the substrate. For now however, we can perform a 
scaling analysis similar to the one we did for the bending 
energy to determine a threshold energy for the pinned- 
to-depinned transition due to in-plane strains, depending 
on the aspect ratio of the perturbation in the substrate. 
In this case the transition is controlled by the ratio of 
pinning energy to elastic energy: 
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where we have used Eqs. (HI) and (11 1. As in the previ- 
ous case, we can argue that the membrane will favor the 
pinned configuration when -jp^ > i^ which gives us the 
condition: 
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(16) 



again consistent with the intuitive picture that shallower 
depressions should favor pinning. The possible equilib- 
rium solutions for the curve of detachment Ix* I are in this 



case given by the extrema of the free energy Eq. (14), as 
in the BD regime, a globally stable solution with \x*\ = 
corresponds to the completely pinned configuration. 
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FIG. 2: Qualitative phase diagram for a membrane on top of 
a patterned substrate of characteristic width L and height 5, 
in units of a the length scale J. The dashed lines correspond 



to the critical lines given by Eq. (171, the solid line is an 
estimated interpolation. Note that this phase diagram is not 
valid for ID geometries, as discussed in the text. 



Equations ( 13 1 and ( 16 1 define two lines of critical val- 



ues given by S^ = S{Lc), which mark the transition from 
pinned-to-depinned in the parameter space of height and 
width of the substrate's profile. Note that while in the 
BD regime the dependence of Sc on the critical width 
Lc is quadratic (see Eq. (13l), in the SD regime this de- 
pendence is linear. In the intermediate region hence it 
is to be expected a crossover between the two critical 
lines. These considerations allow us to construct a qual- 
itative phase diagram. For this it is useful to consider 
the dimensionless quantities S^S/I., L^L/H that give 
the height and width of the substrate's profile in units 
of the length scale 5 of the BD to SD regime crossover 
defined in Eq. (Isl. We can then write the critical lines 
as: 



Sr^hl^Ll 5c <1 



S,= 
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(17) 



Lc 5'c > 1 



The qualitative phase diagram is shown in Figure [2] 

Near the critical line the free energy Eq. (fTl) can be 
written as a Landau functional of the order parameter 
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where A2 is a positive constant (in accordance with 
|xq| = being always a local minimum) and the coef- 
ficients Ai, i = 3 , ...n are functions of the control param- 
eters S, L. As usual, the expansion is cut at order n > 3, 
being An the first non-negative coefficient. The powers 
appearing in the expansion are dictated by the symme- 
try of the system, for cylindrically symmetric geometries 
only even powers are allowed. In the following sections 



Eq. (10) take the form: 




FIG. 3: Membrane on top of a substrate with a depression. 
The figure is axially symmetric with respect to the vertical 
axis through the center of the substrate's depression. R in- 
dicates the radius of detachment. The dashed line represents 
the approximation used in Section |IV| for when the in-plane 
modes are taken into account. 



we will re-obtain these results in an analytical fashion for 
certain simple geometries of the substrate. 



III. DETACHMENT DUE TO OUT-OF-PLANE 
MODES FOR RADIAL SYMMETRY 
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(21) 



Applying the boundary conditions Eq. (21 1 over the gen- 
eral solution Eq. (20 1 we obtain: 
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(22) 



The second equation determines the radius of detach- 
ment, but also imposes a condition over the substrate 
profile for a non-trivial solution to exist (note that there 
is always a solution with R = 0). As we pointed out in 
Section |TT] the radius of detachment R corresponds to ex- 
trema of the total energy of the membrane, which within 
the present approximation consists of 



-^BD (-R) — Epin{R) + Ei^{R) , 



(23) 



with Epin the pinning energy and E^, the bending energy. 
In cylindrical coordinates these are given respectively by 



As we stated in Section |IT] for ID geometries the so- 
lution obtained by only considering the bending rigidity 
term in the elastic free energy, is exact. For three di- 
mensional (3D) geometries this is an approximation that 
works well for small height fluctuations of the substrate. 
Our aim in this section is to obtain analytical results in 
this limit to obtain a qualitative understanding of the 
depinning process. Analytical results can be obtained 
for certain simple geometries, we will restrict our anal- 
ysis to cylindrically symmetric cases. We consider first 
a substrate with an axially symmetric depression s{r) 
as shown in Figure p] The bi-harmonic equation nh in 
cylindrical coordinates, assuming a rotational invariant 
case, is given by: 



-dr + d^] h{r)^0 



(19) 



Epin{R) ~ls I TTrdr , 
Jo 



E^{R) = TTK / rdr 
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dr + dl j h{r) 
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(25) 



In both these expressions the energy is measured from the 
totally pinned configuration. The total energy E'^^{R) 
allows us to determine the stability of the solutions R = 
and Eq. ( 22 1 . This is simply exemplified for the case of 



a parabolic well. 



Parabolic Well 



which has the following general solution: 



hir) = Ho + H.logi 



H- 



2 2 

r 



Hsr^logi 



(20) 



For a depression, if we assume that the membrane de- 
taches from the substrate homogeneously at a circumfer- 
ence of radius R (to be determined) , the solution Eq. ( 20 ) 
is valid for < r < i? and hence it has to be regular 
at the origin. Hi = H^ = 0. The radius r = R gives 
the parametrization of the curve of detachment x* intro- 
duced in Section [IT] The boundary conditions Eq. (Is])- 



In the case of a parabolic profile s{r) = So + 2 J'^, we 



see that (22 1 implies that h(r) = s{r) and there is no 



solution for a partially detached membrane. If we allow 

for a quartic term s{r) = Sq + ^1 

that the detachment radius is given by: 



§jr^ then we obtain 



fif = 



S4 



(26) 



where we have taken the minus sign in the second equa- 
tion of (22) corresponding to the fact that the curva- 



ture of the detached membrane is smaller than the one 



of the substrate. Therefore S4 needs to be a positive 
quantity for a solution to exist. From the total energy 
E]^'^{R) given by Eq. p3| , it is easy to show that the 
solution given by Eq. ( |26p corresponds to a maximum 
of the energy profile and hence it is an unstable equi- 
librium solution while i? = is a metastable minimum. 
£'3°'(i?)— )■ — 00 for R—>-oo and therefore R^ signals the 
energy barrier for total depinning which is always the 
stable configuration. This is however a construction of 
the unbounded quartic profile we have chosen for the 
substrate. In the next sections we will study in detail 
a more physically sensible profile: a Gaussian depression 
or protrusion. 



B. Gaussian Depression 

A more realistic landscape for the substrate is the case 
of a Gaussian depression 



s(r) = G, 1-e 



(27) 



For this geometry, the curvature of the substrate varies 
from positive to negative along the radial coordinate 
and therefore we have to allow for both signs ±-v/^^ in 



Eq. ( |22| . However when the condition is applied to this 
particular shape we obtain: 
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that is, only the positive sign leads to the existence of a 
solution since we have assumed Gg > 0. The solution for 
the membrane's profile then is given by 



h{r) = 




1 



20-2 



Q<r<R 
r> R, 



with R given by ( 28 1 



(29) 



We apply now this solution to the particular case of 
graphene on top of a Si02 substrate. As we did in 



Eq. (17 1, in what follows we will treat all length quanti- 



ties as dimensionless, given in units of the characteristic 
length i. defined in Eq. (Isl. We can consider, as an ex- 
ample to illustrate the solutions given by Eq. (28 1 and 



Eq. ( 29 1 , a particular substrate depression of amplitude 
Gs = 1 and width a = 2. We obtain two possible solu- 
tions for a partially detached configuration: i?i w 0.45, 
and i?2 ~ 2.75. In Figure[4]we depict the graphene mem- 
brane profile solutions that correspond to this particular 
configuration. 

To study the stability of the obtained solutions, we use 
the total energy which, by Eq. (23 I is given byE3 
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FIG. 4: Graphene depinning from a Gaussian depression in 
the BD regime. The plots show the profile of the graphene 
membrane h{r) (Color) and of the substrate s{r) (Black), in 
units of 5, as a function of the distance from the center of the 
depression in units of the depression characteristic width a. 
The figures have axial symmetry. Main figure: stable solution 
with the graphene sheet almost completely detached (-R2 ~ 
5.5). Inset: Unstable solution with very little detachment of 
the graphene sheet (7?i « 0.9). Height profile as a function 
of the radial distance to the center of the depression, r. We 
have taken Gs = 1 and a — 2. All quantities are in units of 
the characteristic length f ~ 1 A. 



The total energy E'^'^{R) corresponding to the substrate 
profile shown in Figure |4] is given in Figure [5] as a func- 
tion of the detachment radius R. We see that the com- 
pletely pinned situation (i?o = 0) is a metastable state 
with a very small energy barrier to overcome to reach 
the true minimum i?2- The solution i?i corresponds to 
a maximum of the energy and hence it is an unstable 
configuration. 

By re-scaling the radius of detachment R by the 
depression width a, R = R/a, the rescaled energy 

-^BD I ^) /'^^ depends only on the ratio Gg/cr^ and hence 
the results can be expressed in an universal manner. 
In Figure 6 we show the re-scaled energy S^^* ( R j /cr^ 

landscape for various values of GgjcP' ■ As expected, the 
global minimum corresponding to the partially detached 
configuration evolves into a metastable state as Gs/cr"^ is 
decreased, and disappears completely for small enough 
Gs/o"^- As discussed in Section Kll the threshold value 
for a stable pinned configuration given by Eq. ( 17 1 is 

Gs/CT^ = 



(30) 



! 0.05. As it can be seen from the figure, 
this estimated threshold is in excellent agreement with 
the exact results. 

In Section |TT we stated that the length of the detach- 
ment curve |x* is the natural order parameter that con- 
trols the pinned-to-depinned phase transition of the sys- 
tem. This can be easily seen now from Figure [6] Given 
the cylindrical geometry of the problem, the length of 
the curve x* is given by |x* | = 2'kRI and hence we can 
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FIG. 5: Main figure: Total energy as a function of the de- 
tachment radius R in units of the characteristic width (t, 
for a Gaussian depression in the bending rigidity dominated 
regime. Inset: Close up showing the metastable solution at 
i?o = and the unstable solution R\ . Results for Gs = 1 and 

(7 = 2. 



take R as our order parameter. As discussed, from Fig- 
ure |6] we see that the minimum at finite K evolves into a 
metastable state that disappears for shallow enough de- 
pressions, while i? = is the true minimum in this case, 
indicating that the transition is a first order one. This 
can be seen in an alternative way by following the evolu- 
tion of the order parameter R. From Eq. ( 29 ) again we 

in 



note that R is controlled solely by the ratio Gs/cr'^, 
agreement with Eq. ( 17 1 and with the universal form of 

the rescaled energy i^^"' ( ^j /ct^. The behavior of .R as 

a function of GsjcP' is shown in Figure m where we see 
that R jumps from i?, = to a finite value at a critical 
value Gs/cr^l ~ 0.2. The figure also shows the spinodal 
point, that is, the value Gs/cr^| ~ 0.05 at which the first 
metastable solution appears, in agreement with the esti- 
mated threshold value. The difference between Gc/cr^l 
and Gg/cr^l shows that the system in the BD regime 
limit is strongly hysteretic. The critical point Gs/cr^l 
and the spinodal point Gs/ct^I can be estimated from 



an expansion of the free energy Eq. ( 30 1 in the re-scaled 
order parameter R: 



T^Tot 




(31) 



This expansion can be identified with the Landau expan- 
sion Eq. ( 18 ) and it assumes that the re-scaled order pa- 



rameter R is small, and hence (assuming a weak first or- 
der transition) it is valid near the critical point Gs/cr^ . 
By minimizing Eq. ( 31 1 it is easy to see that the condition 



for the existence of metastable solutions with i? 7^ is 
given by (Gs/cr^) > (27/4) 7s « 0.12, a value which is 
of the order of magnitude of the spinodal point obtained 



FIG. 6: (Color online) Energy landscape as a function of de- 
tachment radius in the bending rigidity dominated regime 
for a Gaussian depression with varying %^. Solid (blue): 
ff|^ = 0.25; dashed (purple): ff J2 ~ 0.08; dash-dot (or- 
ange): ffig « 0.06; dotted (red): ffj^ Ri 0.03. The energy 
presents a minimum for the detached configuration for cases 
1, 2 and 3 (being this last one metastable), while for case 4 
the energy is a minimum only for the completely pinned con- 
figuration, in agreement with the threshold value discussed in 
the main text. For cases 1, 2 and 3 the completely pinned 
configuration is a local minimum with a low energy barrier, 
not visible due to the large scale of the plot. Note that a and 
Gs are dimensionless: a , Gg ^ tr/C, Gs/J. 
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FIG. 7: (Color online) Order parameter R = R/a as a func- 
tion of the ratio Gs/o-^ of a Gaussian depression in the BD 
regime. The dashed (red) line indicates the spinodal line, 
while the solid (blue) line corresponds to the true transition. 



exactly in Figured For (Gs/a^)^ (> 27/4) 75, the ex- 
trema condition oE^'^/dR = in expression Eq. (31 1 



renders Rq = plus two real positive roots in agreement 
with the energy profiles presented in Figure [6] for the ex- 
act solution Eq. ( 30 ) . 
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FIG. 8: Membrane on top of a substrate with a protrusion. 
The figure is axially symmetric with respect to the vertical 
axis through the center of the substrate's bump. As discussed 
in the text, the membrane is shown as depinning from the top 
of the bump, and re-attaching at a radius L. 



brane follows the substrate for < i? < r and then 
detaches "forever" is not possible. The most general so- 
lution is to assume that there is a radius of detachment 
R and a radius of re-attachment L, with R < L. The 
gain in pinning and bending energies with respect to the 
totally attached configuration in this case are given by: 



(33) 



Epin{R) R 


i js / Trrdr , 
Jr 


= TTK / rdr 
Jr 


(-dr + df)h{r) 
\r J 


— TTK, rdr 
Jr 


Vl o\ 1 

-dr + dl s{r) 
\r J 



(34) 



C. Gaussian bump 

For a Gaussian protrusion 

s{r) = Gse~^^ 



(32) 



the general solution Eq. ( 20 1 holds for r > R. Since 



the origin is avoided, Hi and H3 can be different from 
zero. This solution however diverges for r~^oo unless 
h{r) = const, which in turn cannot satisfy h'{R) = s'{R). 
Hence a kind of solution for which the graphene mem- 



Intuitively it is easy to see that if the bump is very 
pronounced, R has to be approximately zero, otherwise 
the bending energy cost is too high. A sketch of the 
system is shown in Figure [8] Assuming this type of con- 
figuration, the general solution Eq. ( 20 1 holds now for 



< r < L, with Hi — for it to be regular at the origin. 
The constant H^ in this case is allowed to be finite since 
the contact force is acting at r = ^^. Imposing continu- 
ity of the solution and its first derivative at r = L, and 
h{0) — s(0) (note that h{r)—i'Ho for r— ^0), for the region 
< r < L we obtain for the membrane profile: 



h{r)=G,- 



Gs 
If 



la- 



log(L2)--log(L2)+e^log(L2 



G, 

Z2* 



2^ 



--1- 



' log (r^ 



(35) 

while h(r) = s{r) for r > L. The optimal value of L can be obtained numerically, as previously, by imposing 
the discontinuity of the Laplacian of the solution due to the contact force. However since we are interested in the 
qualitative aspect of the solution, it is simpler to analyze directly the energy profile as a function of the re-attachment 
radius L. As it was the case for the Gaussian depression in Subsection niB[ the energy Eq. (36 1 and the re-attachment 
radius L can be rescaled by the width of the bump to show the universal behavior. 



— -^ — = 7s7r — 



KTT jC 

20-4 



L^ + 6i2+2 9-8e"^ 



16 1 



2 1 
I2 



2a4 



(36) 



with L = L/a. The energy profile Eq. (36) is shown in 
Figure |9]for a bump with G^ = 1 and ct = 2, as a function 
of the re-attachment radius L. From the figure, it can be 
seen that the case of total adhesion of the graphene mem- 
brane to the substrate, in the case of a pronounced bump, 
is a metastable state with a very low energy barrier to 
fall into a configuration for which the membrane attaches 
to the substrate after a finite radius L. 

As we mentioned previously, the solution Eq. ( 35 1 is 



valid in principle for pronounced Gaussian protrusions, 
for which Gg/cr^ ^ 1. However it can be shown that this 
is true for any Gaussian bump. This can be seen more 



rigorously by calculating the most general solution for 
which both the depinning and re-attachment radius are 
finite, and finding the minimum of the energy surface. 
The explicit solution for this most general case is rather 
cumbersome and it is given in Appendix |X] here we show 
a plot of the energy surface profile as a function of both 
detachment and re-attachment radius R and L. As it 
can be seen from Figure |10| the complete solution indeed 
shows that the case R = and finite _L is a minimum for 
the case G^/cr^ = 1, and the same can be shown for other 
aspect ratio protrusions. From Eq. (36 1 it is evident that 

There- 
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FIG. 9: Main figure: Total energy as a function of the re- 
attachment radius L for a pronounced Gaussian bump in 
the bending rigidity dominated regime. Inset: Close-up near 
the metastable configuration corresponding to a completely 
pinned membrane. We have taken Gs = 1 and a = 2. 
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FIG. 10: Energy profile as a function of detachment radius 
R — R/a and re-attachment radius L = L/a with Ga/cr^ = 1 
for a Gaussian bump. The global minimum corresponds to 
R = Q and finite L. 



scaled energy profile projection onto the i? = plane, 
E]^"^ (L) / s"^ , for varying Gg/cr^ is shown in Figure [ill 
showing the crossover from the pinned to the partially de- 
tached configuration for increasingly pronounced bumps. 
As in the case of a Gaussian depression analyzed in the 
previous subsection, the case of a finite re-attachment ra- 
dius L is the energy minimum for Gs/<j'^ ^ 0.05, while 
in the opposite limit the minimum corresponds to L — 0, 
that is, for smooth bumps the membrane minimizes its 
energy by conforming completely to the substrate. 

Similar results to those obtained in Subsection IIIIBI 
for the behavior of the re-scaled order parameter i? as a 
function of Gg/cr^ can be obtained here for L, showing a 
first order phase transition between the pinned and de- 



FIG. 11: (Color online) Energy landscape as a function 
of re-attachment radius L = L/a in the BD regime for a 
Gaussian protuberance of varying ratio ^. Solid (blue): 
%^j — 1; dashed (purple): %^L = 0.25; dash-dot (orange): 
^ « 0.l| ; dotted (red): ^ « 0.04J^. The energy presents 
a minimum for the partially detached configuration for case 
1 which disappears completely for case 4, going through a 
metastable state for case 3. Note that a is dimensionless: 
a -^ a/l 



pinned phases. 



IV. DETACHMENT DUE TO IN-PLANE 
MODES FOR RADIAL SYMMETRY 



In the previous sections we studied the detachment of 
a graphene membrane from a patterned substrate due 
to the bending rigidity term in the free energy Eq. (pi. 
This approximation is widely used, but, as we showed, it 
is valid for relatively small fluctuations of the substrate 
landscape for a 3D pattern. In this section we consider 
the less studied case of depinning due to in-plane modes, 
for which the free energy Eq. (fTl) is approximated by 
Eq. (14 1. As before, we will restrict our study to cases 



that allow for an analytic solution, in particular a sub- 
strate with radial symmetry. In cylindrical coordinates 
this elastic energy is given by: 



Eel = - I 2nrdr 



ji / 2TTrdr 



drUj. H — - -\ — {drh) 
r 2 



drUr + - {drh) 



-I 2 



ji I 2TTrdr I — ) 



(37) 



where Ur is the radial component of the in-plane displace- 



ments, and Epin was defined in (24). In the following 



subsections we will obtain results for a Gaussian depres- 
sion and protrusion. 
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A. Gaussian Depression 



differential equation to solve is given by: 



In this section we consider again a Gaussian depression 
given by Eq. ( 27 1 , over which there is a membrane par- 



tially attached. For radius greater than a radius R the 
membrane is pinned to the substrate and follows its pro- 
file while for < r < R the membrane is completely de- 
tached. When analyzing the effect of the in-plane modes 
we encounter an added complication which is that we do 
not know the height profile of the membrane for the de- 
tached region, since this would imply to solve the problem 
completely by treating the full coupled non-linear differ- 
ential equations in both h and u fields resulting from 
minimizing Eq. M . Here we consider as a first approxi- 
mation that the graphene membrane remains fiat within 
the detached region as shown in Figure [3] Hence the 



- (A + 2^) df.Ur + 







dr 



0<r<R 



(38) 



(A -I- 2n)drs {dfs) + t^ (drsY r>R 

Our ansatz corresponds to a membrane profile given by: 



h{r)^ 



ho 

s{r) 



0<r<R 
r>R 



(39) 



with s{r) given by Eq. (27). The general solution of 
Eq. (38) is given by (see App. [B|: 



Ur{r) 



4(t2 (A+2p)' 
2 r 



Gt 



(A+p) 1 
(A+2^) r 



1 Gj (A+At) 
r 4ct2 {\+2fj.) 



(i?2 



0<r<R 
r>R. 



(40) 



The radius R of detachment can be found by finding the extrema of the total energy _EJ°'(i?) = Epin{R) + Eei{R) 



where Epin is given by Eq. (24) and we measure the elastic energy Eq. (37) from the totally attached configuration 



The total energy as a function of detachment radius can be calculated to be: 



E^^^'iR) 



2i?2 



0-4 8(A + 2fi) 



M^ 



(41) 



Minimizing £'J°'(_R) renders a solution with Rq = 
which is a local minimum and corresponds to the mem- 
brane completely pinned, and the following transcenden- 
tal equation for the equilibrium detachment radius: 



Gt i?2 , (A + ^) _2if 



2ct4 fj2 



ix + 2^iy 



(42) 



Again, as in the BD dominated regime case of Section pII| 
we see that by re-scaling both the detachment radius and 
the total energy by an overall factor given by the depres- 
sion width a, R = R/a and £^J°*(i?)/cr^, the re-scaled 
energy shows universality. In this case, and in agreement 
with Eq. (17), the system is controlled by the ratio Gs/u, 
opposed to the dependence on Gg/cr'^ found for the BD 
regime. 

We can apply our results to a graphene membrane on 
top of a Si02 substrate as we did in Section |III[ Taking 
the accepted values for room temperature for the Lame 
coefficients of graphene, /x w 10 eV A~^, A ~ 2 eV A~^ 
ESlj and G^ — 5 (in accordance with the validity of our 
approximation) and cr = 4 (note that we are still work- 
ing in in units of the scaling length J) we get two possible 
depinning radius, i?i « 0.1 corresponding to an unstable 
minimally detached configuration, and i?2 ~ 8.4 which 
is the stable, global minimum solution. The total en- 



ergy i?J°' is plotted as a function of detachment radius 
showing the different equilibrium 
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R = R/a in Figure 
solutions. The energy barrier to be overcome for detach- 
ment from the metastable equilibrium configuration at 
Ro = is very small as can be seen from the inset in Fig- 
ure [12] Note also that the minimum at R2 is very shallow 
in comparison to the energy scale, as shown in the inset 
of Figure [12] and hence any fiuctuation could lead to the 
graphene membrane to be detached at a radius R > R2, 
and therefore closer to the fiat configuration. This ef- 
fect is less pronounced as the width of the depression is 
increased. 



In Sectionjil] Eq. (17), for the SD regime we predicted 
a critical value for the ratio Gs/a above which the stable 
configuration is that of the membrane partially detached 
from the substrate. Taken the values for graphene dis- 

/ \i/4 

cussed above, this ratio is given by (4j^J w 0.1. 

This estimate is in good agreement with the exact re- 
sults for the re-scaled energy profile i?J°'(i?)/(T^ depicted 
in Figure [13] as a function of the re-scaled detachment 
radius R and varying Gg/cr. 



As discussed in Section |IIH the radius of detachment 
R can be taken as the order parameter of the problem. 
We show the behavior of i? = R/a as a function of Gg/a 
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FIG. 12: Energy profile for a Gaussian depression as a func- 
tion of detachment radius K = R/a for Gs = 5 and cr = 4 
as discussed in the main text. The insets show the local min- 
imum that corresponds to the totally pinned configuration, 
and the global minimum corresponding to the partially de- 
tached membrane. 



in Figure 14 Again, a pronounced jump in i? to a finite 
value with increasing Gs/<J is observed, at a critical value 
Gs/cr|^ in agreement with the predicted threshold for the 
transition. The spinodal line, also shown in Figure [141 is 
basically indistinguishable from the true transition and 
hence in the SD regime there is almost no hysteresis. 
The expansion of the free energy Eq. (41 1 in powers of 
the order parameter R: 



77'Tot 



'u {r) 



isR^ 



G UiX + lj) 

Cr4 (A-f2/l) 



R^ 

4 



3 



(43) 

gives a good qualitative description of the first order tran- 
sition obtained exactly in Figure [14] Moreover, in this 
case the spinodal point given by the expansion Eq. ( 43 ) 
is given by {GJaf = (44/15) 7s, giving GJa\^ « 0.28, 
in excellent numerical agreement with the exact value. 



FIG. 13: (Color online) Energy landscape in the SD regime as 
a function of detachment radius R — R/a, for a Gaussian de- 
pression of varying ratio ^ . Solid (blue): ^| = 0.5; dashed 
(purple): ^j^ ^ 0.33); dash-dot (orange): ^j^ = 0.25; dot- 
ted (red): %^L ~ 0.17. The energy presents a minimum for 
the detached configuration for cases 1, 2 and 3 (being this 
last one metastable), while for case 4 the energy is a mini- 
mum only for the completely pinned configuration, in agree- 
ment with the threshold value discussed in the main text. For 
cases 1 and 2 the completely pinned configuration is a local 
minimum with a low energy barrier, not visible due to the 
large scale of the plot. The region of metastability for the 
depinned configuration is very small, as seen in this plot and 



also in Figure 14 Note that a is dimensionless: a — > a/t 



general solution valid for small protrusions (BD regime) , 
Eq. (20 1 ■^'*. The details for the solution using this ansatz 

Imposing the boundary con- 
leading to the 
5 and 



are given in Appendix p] 

ditions results in two possible solutions, 

two energy profiles shown in Figure [15] for Gg 

CT = 4. Although these solutions differ for the metastable 

or unstable regions (a construction of the approximation 

involved), they coincide for the minimum and hence the 

stable re-attachment radius L is uniquely defined. 



B. Gaussian bump 

Following similar manipulations to the previous sec- 
tion, we can calculate the solution for the radial com- 
ponent for the in-plane displacements in the graphene 
membrane due to a Gaussian protrusion parametrized by 
Eq. 



(32 



section 



Given our findings for the BD regime in Sub- 



III C| we consider a configuration of membrane 



on top of the substrate that is pinned at the very top 
(detachment radius R — 0) and re-attaches at a radius L 
as shown in Figure [8] This takes into account the ener- 
getic cost of bending. As in Subsection IIIB we have to 



make a sensible approximation for the unknown profile 
of the membrane on the detached section. We hence ap- 
proximate the detached profile of the membrane by the 



V. DISCUSSION 

In this work we have analyzed the possibility of depin- 
ning of a membrane on top of a patterned substrate. We 
have studied simple configurations of the substrate that 
allow for analytical solutions, in the two relevant limits 
for the problem: the bending rigidity dominated regime, 
valid for small corrugations of the substrate, and the 
elastic, strain dominated regime which holds for larger 
corrugations. We have shown that in both limits, the 
energy cost of either bending or stretching can cause the 
membrane to prefer to detach from the substrate. The 
particular results confirm the more general picture we 
sketched in Section [ll] in which we obtained a qualitative 
phase diagram for the system presented in Figure [2] This 
phase diagram presents two phases, a completely pinned 
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FIG. 14: (Color online) Order parameter i? as a function of 
height Ga of the Gaussian depression in units of a in the SD 
regime. The dashed (red) line indicates the spinodal line, 
while the solid (blue) line corresponds to the true transition. 
In the SD regime these two lines are indistinguishable. 
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FIG. 15: (Color online) The two possible solutions for the 
energy profile as a function of re-attachment radius L = L/a 
in the SD regime, as calculated in AppendixlB] We have taken 
Gs = 5 and cr = 4. 



phase in which the membrane follows the profile of the 
substrate, and a detached phase in which the membrane 
prefers to depin partially at some optimal detachment 
curve X*. We have shown that an order parameter can 
be constructed from this detachment curve which allows 
us to write the problem in terms of a Landau energy func- 
tional Eq. (18). By scaling arguments we were also able 



to obtain the critical lines ( 17 1 in both BD and SD limit 



ing regimes and found excellent agreement with analyti- 
cal calculations. The critical line in Figure [2] represents 
a first order phase transition, as argued in Section |TT] and 
shown explicitly for specific geometries of the substrate. 
We showed that the energy of the system shows univer- 
sality, and depends on the ratio of S/L or S/L"^ for the 
SD and BD regimes, where S is the height of the sub- 



strate's corrugation and L its characteristic width. We 
obtained the critical and spinodal points exactally for the 
analytically solvable cases, and showed that the Landau 
energy functional gives a very good estimate to locate 
these points. 

The depinning process is dependent on the aspect ratio 
of the spatial perturbations of the substrate, the elastic 
parameters of the membrane, and its interaction with 
the substrate. The interaction between graphene and 
different types of substrates is not well known. Order 
of magnitude estimates for different mechanisms '^ sug- 
gest that the interaction coupling 75 ~ 10~^ — 2 meV A~^ 
(note that the interaction between two graphene layers in 
graphite is 20-30 meV A"^'^^. Within this work we have 
set 75 = 2meV A~^, its most conservative value, and 
hence the obtained values for depinned configurations are 
underestimations. For this value of 75 we find that, in 
the elastic regime, the depinning of graphene becomes 
relevant for height corrugations such that S/L > 1/10. 
In the bending regime, the condition depends on the to- 
tal area of the corrugation: S/L'^ > 1/20 A . These 
obtained values for possible depinning are comparable to 
measured corrugati ons in free standing graphene '^, and 
in graphene on Si02'^^^. Hence, regions where graphene 
is detached from the substrate may be found in samples 
on Si02, in agreement with the observations reported 
irp3. Although these conditions were obtained for the 
two limiting regimes, in real life both effects are present. 
In general, corrugations of all scales are ubiquitous due 
to the intrinsic roughness of the substrate, and both the 
bending rigidity and the in-plane strains of the graphene 
membrane will contribute to its depinning. Random con- 
figurations of the substrate within a mean field model as 
presented in Section |TT] could be treated by adding noise 
to the system, considering the parameters in the Landau 
free energy Eq. ( 18 ) as random, in the spirit of random 
mass theories. 

The presence of corrugations in graphene, either intrin- 
sic or substrate-induced, can lead to diverse experimental 
consequences. Corrugations in graphene are associated 
with gauge fields that couple to the Dirac electronJiSl. 
These gauge fields generate an effective magnetic field 
that can affect the transport properties of graphene. We 
have shown however that, for graphene, the scale of the 
corrugations for which the bending rigidity is relevant 
is rather small. As we saw in Section |III| the BD to 
SD regime crossover length for graphene is f sa 1 Aand 
hence depinning due to the strain energy cost is to be 
expected. The effective magnetic field is related to the 
strain through the relation 



B 



f3a 
a L 



(44) 



where 4>q « 10^^^ Wb is the quantum of magnetic flux, 
a = 1.42 A is the lattice parameter and /3 « 2 gives 
the change in the hopping parameter between nearest 
sites for a Dirac electron due to the deformation of the 
lattice. The corrugation of the membrane determines 
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L, the characteristic width of the corrugation, and the 
strain a, which for simplicity we assume to be con- 
stant. We can roughly estimate the maximum magnetic 
field that the graphene membrane can experience do to 
the corrugation of the substrate. The strain of a cor- 
rugation of height s and width L scales as s^ / L? and 
hence, from Eq. (16 I, the maximum strain that graphene 
can support is am ~ 2%. As we showed, beyond this 
point the membrane relaxes by depinning partially from 
the substrate and lowering the strain. Taking a phys- 
ically relevant corrugation width of L « 100 nm, a^ 
corresponds to a maximum effective magnetic field of 
Bjn « 2 T. This order of magnitude indicates that, in- 
deed, the magnetic field due to the induced corrugations 
can have a sizeable effect in the transport properties of 
graphene ^ '^^ . Our results indicate that rougher sub- 
strates could in fact lead to flatter configurations of the 
graphene membrane after annealing, due to the impos- 
sibility of graphene to conform to pronounced depres- 
sions or bumps. This could result, in a counterintu- 
itive fashion, in greater mobilities for graphene on top of 
very rough substr ates, du e to a decrease in impurity and 
phonon scattering^SESEIl On the other hand, for these 
kind of substrates the graphene membrane would be al- 
most suspended and therefore prone to the excitation of 
flexural modes which contribute to the resistivity^^'^^. A 
more unexplored path is the possibility of controlling the 
pinning or not of graphene to the substrate by tuning 
the different metastable states, which could be realized 
by modulating the gate electric field or by applying exter- 
nal pressure. As we showed, the system can present hys- 
teresis and irreversibility. Our results can be also help- 



ful to the understanding of the ubiquitous formation of 
graphene bubbles on many types of substrates'^*^ ''■^ . 

To conclude, we list the limitations of our model. Our 
work is based on a continuum approach and hence it 
breaks down for lengths of the order of the lattice spac- 
ing. However amplitudes of the order of the lattice spac- 
ing can still be well described by the continuum model 
as shown in Ref. [251 As we pointed out, the phase di- 
agram presented in Figure [2] is valid for 3D profiles of 
the substrate. The crossover to ID geometries would in 
principle imply the disappearance of the linear critical 
line valid for the SD regime, since the BD results are 
exact for ID. The same holds for exact calculations of 
Section IIIII and Section 11 V I which have been done for a 



isotropic perturbations. The effect of the lack of radial 
symmetry remains to be explored. Lastly, it is possible 
that the interaction between the graphene layer and the 
substrate is not uniform, due to the presence of charges 
and other defects within the substrate. The modeling of 
this kind of potential goes beyond the scope of this paper. 
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Appendix A: Full solution for a Gaussian bump: flexural modes 



In the main text we have analyzed the stability of a membrane on top of a Gaussian bump, which is pinned at 
the very top of the bump an re-attaches to the substrate at a radius L. This was an assumption based on energetic 
arguments. The most general solution can have the membrane conforming to the substrate from the top up to a 
finite depinning radius R, and then re-attaching at a radius L. In this case the general solution Eq. (20 1 is valid for 
R < r < L and all the coefficients Hi, i = 0, ...3 can be finite. These can be found by imposing the continuity of the 
solution and its first derivative both at R and L. Defining 



n 



G..e 



(L2 _ R2f _ 4X2^2 (log I _ log ^Y 
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we have: 



{L^ - R^) (e^L^ - c^bA a"^ - 2e^L^R^ {R^ + 2a^) \ogL' 



HL^R^ log R 2e^cr^ + e^ {L^ - R^ - 2a^) - 2e^L^R'^ (L^ + 2a^) logR 



+ nL^R^ log L 
Hi = -HL^R^ ff 



H. 



2'HL^R^ 



'2H 



2e-^a^ ^e"^ {-L^ + R^ -2a^) +2{e^ {L^ + 2cr^) + e^ {R^ + 2<t^) ] logR 
:^ - e^] (L^ - i?2) + (2e^ (L^ + 2a'^) - 2e^ {R^ + 2a^)] logL 
^ {L^ + 2a^) - e^ (i?^ + 2a^)] logR 
(e& -e^\ {L^ - R^)a'^ ~2e^L^R^ log L^ 



2HR^logR 



2e^^(j^ + e^^ (L^ - R^ -2a^) -2e^^L^logR 



21-LL^logL 2e^a'^ + e^ i-L^ + R^ - 2a'^) +2{e^ + e^ ] R^logR 



H3 = -H {L^ - R^) et^ (l2 j^ 2(7 



(i?2 + 2o-2) 



2HL^R^ 



(6 2^2 — e 2^2 j log L — I e 2<t2 — e 2<t2 j log R 



(Al) 



The total bending energy can be calculated as in the main text taking into account that we now have three different 
regions of integration. The result, measured from the totally pinned configuration, is given by: 



E^ = KTrV. 
+ kttT-L 



{L^ ~ R^ 
(L^^R^) n 






2ff4 

mTn ^^\'!^^K ^ [R^^ + 2a^{3R^ + 9R^a^ + 8a^) - L^{R* - 2R^a^ + 2a4)] 



in-Ky. 



2KTTn 



2kt:'H 



2^4 
L2i?2(logL-logi?) L-^+R 



(A2) 



e^:;^(L2 + i?^+4cr^) 



L2i?2(logL-logi?) 4 



[4cr2(i2 j^ 2cr2) + {L^ + 22.2 0-2 + 2a^) (logL - logi?)] 



L^ R^ {log L- logR) 4 



e^ [-4cr2(i?2 _^ 20-2) + (i?4 _^ 2i?2cr2 + 2(7^^) (log i - log i?)] 



This bending energy together with the contact energy cost gives the plot shown in Figure [TO} 



Appendix B: Solution for a Gaussian depression and bump: In-plane modes 

1. Gaussian depression 



We will call u^: and u^ the solutions of ( 38 1 for < ?■ < i? and r > R respectively. For < r < _R ( 38 1 is 



homogeneous and has a general solution of the kind: 



Ao 



Air . 



The constant Aq — for the solution to be regular at the origin, and Ai is determined by the boundary conditions. 
For r > R the homogeneous solution is given by 



Co 



Cir. 



but in this case we get Ci = by imposing that the displacements are at infinity. Cq is determined by boundary 
conditions. To obtain the general solution of (38 1 for r > i? we have to add a particular solution. This can be obtained 
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by the ansatz solution: 



= ^a„ 



Substituting (|Bl| into the second hne of (|38| and using (|39|) we obtain 

4 



J2a^ 



„in+2 



;{m+ l)r™ + (m^^ - l)r 



m-2 






From here we see that there is a possible solution with m = ±1. By substituting in (B2| we obtain: 

G2 (A + m) 



fli 



4cr2 



fl-l 



4 (A + 2^) 



(Bl) 



(B2) 



(B3) 



The coefficients Ai and Co are determined by imposing the continuity of the solution u^ (R) = u^ {R) and of the 
in-plane stresses at R: 



Grr — A 



drUr H \ (drh) 

r 2 



2^ 



drUr + - {drh) 



(B4) 



and hence the full solution is given by ( 40 ) in the main text 



2. Gaussian bump 

As discussed in the main text, for a bump we approximate the membrane's profile in the detached region by the 

Ho 



general solution given in the BD regime Eq. ( 20 ) 



h{r)^Gs+^-fr'' + H^r''\ogr 



(B5) 



where we have set Hq = Gg and _ff i = as explained in Subsection |IIIC[ If we denote the re-attachment radius asL, 
then the in-plane radial displacement w,. is given by: 



- fA + 9n\ (rl^. + ^'■"'- _''r\_{{\ + 2^i)^rh (9» + ^ {Orhf 0<r<L 

[^ + 2^^)^d,Ur+ ^ r^)-{(X + 2^^)^rs{^'rs) + >^{^rsf r>L 



(B6) 



together with the appropriate boundary conditions. The solutions for < r < L and r > L are given respectively by: 



u<{r) =Cir 



+ 



r^ [-2HiH2{\ + m) + 2gf (A + 3/^) + Hj{\ + ^)] 
16(A + 2/i) 
'H2 log (r) [H2{\ + ^i)-2{Hi+ H2 log (r)) (A + 3/i)] 



u>{r) 



Co Gie 



2„-i2 



r 4(A + 2^)r 



4(A + 2^) 
(A + /i) + (A + 2n) 



(B7) 



The parameters Hi, H2, Cq and Ci are fixed by the boundary conditions, continuity of the solution Eq. (B7l and its 



first derivative at r = L and continuity of the in-plane stress Eq. (B4|, plus continuity of the fiexural field h{L) ~ s{L) 
Imposing these result in two possible sets of solutions: 



Co = 

( 
Ci=- 

Hi 
Ho 



Gle-^{X + /i) L* + AL^cr^ + 8cr'* - Se^cr" + 4e^cr 



16(A + 2/i)CT4 

'' 4L2(A + 2/i)cr4 

2G,e"^ (-cr2 + e^cr2 _ ^2 log (^) _ 2o-2 iQg (j;^) _^ 2ei^a2 log (L)) 



(B8) 



£2^2 



( 



G^e 57^ (-L2 -2ct2 +2e5^cr^ 



L2a2 
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Co 
Ci 
Hi 
H2 



'' 16(A + 2^)cr4 

' 4L2(A + 2^)cr4 

2G,e-^ f-cr2 + el^cr^ _^ 2.2 log (L) - 2^2 log (L) + 2e^cr2 ^^g (^) 



(B9) 



L2c 



G^e^^^ ( L2 - 2cr2 + 2e^a' 



The total energy can be calculated from the general expression Eq. (37 1 by use of Eq. (B7) and Eq. (B8l-Eq. (B9l 
resulting in the two energy profiles plotted in Figure [TF 
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Previous works have studied this problem for graphene re- 
stricted to the bending rigidity dominated regime and for 
some particular patterns of the substrate "^tPo i j^ compre- 
hensive treatment of the problem of adhesion of a mem- 
brane to a substrate, in the BD regime, can be found in 
Ref. lT6l37l 

The value of 7s ~ 2meV A~^ is actually an upper limit 
for 75, attained for when there is a layer of water between 
graphene and the substrate-^. Also, repulsive forces not 
considered in Ref. 17|jwould effectively reduce the value of 

is- 

Note that although the total Gaussian curvature of a nearly 
flat 2D membrane is zero, the quartic interaction term gen- 
erated by integrating out the in-plane phonon modes can 



be seen as a long range interaction between the local Gaus- 
sian curvatures at different points in the membrane (see 
Ref. lisp. This long range interaction is the responsible for 
the stability of a low temperature flat phase in 2D mem- 
branes (see Refs. [TI20)| . 

'^^ Compare with the temperature dependent crossover length 
scale obtained by the renormalization of the bending rigid- 
ity due to the in-plane modes, see Refs. I1I38I 

*^ Note that in this general context, the curve x* can be 
disconnex. 

** Note that we have to retain the non-linear term in the 
strain tensor since we are in the limit of large out-of-plane 
fluctuations. 

'^^ For the figures we have used the exact expression for the 



pmnmg energy. 



E„ 



= —jTv J rdr\ 1 + I ^re 



